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Abstract 

In this paper we extend the umbral calculus, developed to deal with difference 
equations on uniform lattices, to g-difference equations. We show that many of the 
properties considered for shift invariant difference operators satisfying the umbral 
calculus can be implemented to the case of the (/-difference operators. This q- 
umbral calculus can be used to provide solutions to linear (/-difference equations and 
^-differential delay equations. To illustrate the method, we will apply the obtained 
results to the construction of symmetry solutions for the (/-heat equation and to 
solve a linear ordinary second order (/-difference equation. 



1 Introduction 



Q-functions appear in many physical problems. They enter in the study of exactly solvable 
models in statistical mechanics [I], in conformal field theory [2|, and are thus very rele- 
vant for applications. For example, g-exponential distributions can be obtained following 
Gibbs' procedure from the stationary conditions on a certain generalized entropy Stan- 
dard g-exponential functions are also used to extrapolate between Fermi-Dirac (q = oo) 
and Bose-Einstein (q = 0) statistics, passing through Maxwell-Boltzmann (q = 1) statis- 
tics m. 

In the case of difference equations one had proved p] that there exists a very powerful 
method for systematically discretizing linear differential equations while preserving their 
properties. Here we extend that method to the case of g-difference equations. We can 
show that many of the properties considered in [3] for shift invariant difference operators 
satisfying the umbral calculus |6p7p8|9pin|ITT] can be extended to the case of the g-difference 
operators considered in Ref . • For any g-difference operator this g-umbral calculus 

can be applied to provide solutions to linear g-difference equations and g-differential delay 
equations. As an illustration, we will apply the method in the construction of symmetry 
solutions for the g-heat equation and to a linear ordinary g-difference equation. 

The paper is organized as follows. In Section 2 we define a g-difference equation 
in p independent variables and, for the sake of simplicity, just one dependent variable 
and characterize the symmetry transformations which will leave the equation invariant. 
Section 3 is devoted to the study of the properties of g-calculus, showing the differences 
and the similarities between differential and g-difference calculus. In particular in Section 
4 we will discuss from an analytic and numerical point of view the simplest g-functions 
which will be used later, in Section 5, along with a few examples. In Section 5 we present 
the symmetries of a g-heat equation and solve a differential g-delay equation using the 
correspondence between g-calculus and differential calculus. Section 6 is dedicated to few 
conclusive remarks. 

2 (/-difference equations and its Lie symmetries 

Let us consider a linear g-difference equation, involving, for notational simplicity, only 
one scalar function u(x) of p independent variables x = (x±,X2, ■ ■ ■ ,x p ) evaluated at a 
finite number of points on a lattice. Symbolically we write 

E N (x, T a u{x), T a ^A Xi u(x), T ai ^A Xii A Xi2 u(x),..., (2.1) 
T a ili2 ... iN A Xii A Xh . . . A x . n u{x)) =0, a = (a u a 2 , . . . , a p ), 

where a, aj_i, a ili2 , . . . are multi-indices, En is some given function of its arguments, 
ii, %2i ■ ■ ■ , in take values between 1 and p. We use the shortening notation T a u(x) = 
{T£*T£* ■ ■ ■ Txp u(x)}, where cij, with i = 1,2, ... ,p, takes values between rrii and rij, with 
rrii, rii fixed integers (rrii < rii), and the individual g-shift operator T£j is given by 

T£u(x) = u(xx,x 2 , ■ • . ,Xi- 1 ,q1 1 Xi,x i+1 , . . .,x p ). (2.2) 
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The other g-shift operators T ai ^, . . . are defined in a similar way. The operator A Xi is 
a gj-difference operator which in the continuous limit when g, — > 1 goes into the partial 
derivative with respect to the x\ variable (few examples of it are given in Section 3). By 
qi we denote a positive parameter which defines a non uniform lattice of the variable 
Xi (i = l,...,p). 

To study the symmetries of eq. ()2.1|) we will use the approach introduced in Ref. 
based on the formalism of evolutionary vector fields for differential equations [T^j. As in 
the case of differential equations the symmetry group of a discrete equation is characterized 
by those transformations of the equation that carry solutions u(x) into solutions u(x). 
Moreover we look only for those symmetries which in the continuous limit go over to Lie 
point symmetries. In such a case the infinitesimal symmetry generators of the symmetry 
group of equation (|2.1)) in evolutionary form have the general expression 



X e = Q(x,u)d u = Ij^U^T^iq^T'A^u-^T^iq^) ) d u , (2.3) 



with £i(x, T a u, {qj} P j =1 ) and 0(x, T c u, {g./}j =1 ) such that in the continuous limit go over 
to £,i(x,u) and (p(x,u), the infinitesimal generators of the corresponding Lie point sym- 
metries. The group transformations are obtained by integrating the differential equation 

du(x) = q(% u(zx\ u(x,g = 0) = u(x), (2.4) 
dg 

where g is the group parameter. Eq. (j2.4j) can be integrated on the characteristics in 
the continuous limit, however in the discrete case the integration of the corresponding 
differential g-difference equation can almost never be carried out. 

Eq. (|2.1j) is of order N in the difference operators and of order i\ + %i + . . . + i n in the 
shift operators. Hence, constructing the following prolongation of X e 

P rX e = J2 Ta QdT«u + J2 Tan Q Xnd T^A Xn u {x) + J^^ 2 ^ 29 ^^^ + ■■■ 

a ii 11,12 

(2.5) 

we find that the invariance condition 

prX e E N \ EN=0 = (2.6) 

must be satisfied. The summations in ()2.5|) are over all the sites present in (12.1}) . By 

Q Xi i, Q Xi i Xi 2, ... we denote the total variations of Q, i.e., 

Q^=A T Xn Q, Q x ^=A T Xn A T Xi2 Q, 

where, in the simple case of the g-right derivative presented below in eq. (j3.1j) . the 
partial variation A Xti is defined by 

^Xi 1 f (^1 j ■ ■ • j ^Cpi u{x\ , . . . , Xp) , A X j u{x\ , . . . , Xp) , . . .) 

= ( 9i —Tjiiii Ifi^i' ■ • • 9*1*^*1) • • • > ^"pi ^(^lj • • • i •Ep)i ^XjU{x\, . . . , Xp), . . .) 
/(*£i> • • • ) u\Xij ■ ■ ■ , Xp) , A X jU\X\i . . . , Xp) , . . .)] , 
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and the total variation A£ by 
A xiJ( x ^ ■ ■ -,x p ,u(xi, . . . ,x p ), A x .u(x u ...,x, 



5 -^PJ 1 



(q 4 —l) Xil [/("^l ' ' ' ' Qii X ii ' " " " ' u(xi, • • • Qii%ii 1 ■ ■ ■ > *£p); ^ay^O^lj ■ ■ ■ Qii%ii 1 ■ ■ ■ 1 -^p)i • • •) 
y ) • • • j "^p> , . . . , j li (xi , . . . , Xp) , . . .)] . 

Notice that expressions (|2.3H2.5j) are analogous to those of the continuous case [T3] and 
can be derived in a similar way [Tfj . 

The group generated by the prolongations also transforms solutions into solutions, 
and A Xi u, A x .A x .u, . . . (up to order N) into the variations of u with respect to the corre- 
sponding Xi. 

The symmetries of equation (|2.1|) are given by condition ()2.6|) , which give rise to a set 
of determining equations for £j and obtained as coefficients of the linearly independent 
expressions in the discrete derivatives A Xi u, A XiX .u,). 

The Lie commutators of the vector fields X e are obtained by commuting their first 
prolongations and projecting them onto the symmetry algebra Q, i.e., 

[X n ,X a ] = [pr 1 X e „pr 1 X a }\ g (2.7) 

V ou ou Ou Xi ou x J 

where the du Xi terms disappear after projection onto Q. 

The formalism presented above may become quite involved, but the situation is simpler 
for linear equations where we can use a reduced Ansatz. In this case we can assume that 
the evolutionary vectors ()2.3|) have the form 



H u - 4>(x,T a ,qj)u J d u . (2.8) 

The vector fields X e can be written as X e = (Xu)d u with 

X = T ^ Qj)A X i ~ <f>(x, T\ q 3 ). (2.9) 



Notice that X span a subalgebra of the whole Lie symmetry algebra (see Ref. 

If the system is nonlinear the simplification (|2.8|) is too restrictive and is almost im- 
possible to get a non trivial result since the number of terms to consider is a priori infinite. 



3 q— Calculus 

In this section we present the generalities of g-calculus jl6j. We will restrict ourself for the 
sake of simplicity to one independent variable. Moreover, in the following we will consider 
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just the simplest g-derivatives (at the right, at the left and symmetric, respectively) 

A* + = -^(?W), (3.1) 

^ = 4-(i-r- 1 ), (3.2) 

Qx x 

Ai = -j-(T„-77 1 ), (3.3) 

where g x . is a real dilation positive parameter associated to the variable x, q l x , i — ±, s are 
given by 

Qx = Qx - 1 q~ = l-— q s x = q x -—, (3.4) 

g^ 

and T x is a g-dilation 

T x f(x) = f(q x x) T~ 1 f{x) = f(x/q x ). (3.5) 
The operator T x is the one-dimensional reduction of eq. (|2.2|) . Formally we have 

r x = e {fc - 1)xfte . (3.6) 

When we do not specify which g-derivative we are using we will write just A x . 

It is easy to see that, due to the form of the g-derivative considered (|3.1I - 13. 3|) . 
the shift operator and the g-derivatives do not commute. So the g-derivative operators 
are not shift invariant operators and do not satisfy one of the basic conditions in the 
umbral calculus 00IE1II3. We can, however, carry out an umbral calculus even if our 
g-derivative are not shift invariant that we will call q -umbral calculus. The g-umbral 
calculus is defined in the same way as the standard umbral calculus, but we do not 
require the commutativity of the delta operators with the shift operators. The absence of 
this property has no consequence in all the results presented in the following. 

We can easily find that 

f Tx for A+ 

[A x ,x] = j T- 1 for A- (3.7) 

I ifeC&r. + T- 1 ) for A* 

If instead of the standard commutator, we consider the g-commutator defined as 
[A, B] q+ = AB - qAB and [A, B] q - = AB - (l/q)AB, we have 

[A x ,x} q = l. (3.8) 

The result (|3.8|) is not valid in the case of the symmetric g-derivative. Moreover, the ex- 
pression ()3.8|) does not satisfy the Leibniz rule and, thus, the expression [A x , x n ] q becomes 
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more and more complicate when we consider n = 2, 3, • • • . So, here in the following we 
will consider just the standard commutator. 

We can define an operator /3 X , depending on T x , in the spirit of umbral calculus [HES], 
such that 

[P x ,T x } = 0, [A x ,/3 x x] = l. (3.9) 

For the three g-derivatives introduced above (|3.1jl we can find the following explicit ex- 
pressions of (3 X : 

/?+ = (q x - l)xd x (T x - l)" 1 = q+xd x (T x - l)" 1 , 

13- = (1 - -)xd.(l - T- 1 )- 1 = q* xd x (l - T- 1 )- 1 , (3.10) 

q x 

PI = (& - ~)xd x (T x - T- 1 )- 1 = q s x xd x (T x - T- l )-\ 
Qx 

These may not be the only possible definitions. 
It is easy to prove that always 

f3 x xA x =xd x . (3.11) 

Moreover, due to the presence of the d x operator in the definition of f3 x , the g-umbral 
correspondence of an explicit x-dependent differential equation will give rise to a q- 
differential delay equation [TTj . 

We can reexpress the functions j3 x as an infinite series in terms of the shift operators, 
thus proving that they commute with the shift operators. From ()3.6|) we get 



xd x = \nT x = ln(l + (T x - 1)) = J2(-l) n+1 — -■ (3-12) 

Qx ~ 1 Qx ~ 1 Qx ~ 1 ^ n 

Consequently 



n+ (T — 1 V™ 

& + = ^£(-ir^— (3 - 13) 

Qx-\^ n + 1 

Similarly from T" 1 = exp (—(q x — l)xd x ) 

xd x = InT- 1 = ln(l + (T- 1 - 1)) = £(-1)"^ >- (3.14) 

Qx - 1 Qx - 1 Qx-1~ n 

we get 



= E(-i) w (rirl " 1)W = E (1 ~I* 1)n - (3-15) 



n=0 ^ n=0 

Finally, for the symmetric derivative, as 



Tx J x = sinh((& - l)xd x ) (3.16) 



we find that 



-| rp rp— 1 i °° (rp rp— l\2n— 1 

xd, = ^-sinh- 1 ^-^ = ^ T ^(-ir 1 a 1 ^ 2 t X J , (3-17) 



where the coefficients C n are given by 

01 = 1 - c - (CTn£fe) ' Vn - 2 ' (318) 

So, we have: 

„s 00 / y _ rp-l\2n 

For all functions / and g, entire in (3 x x, the Leibniz rule takes the form 

[A x ,fg} = [A x J}g + f[A x ,g]. (3.20) 

Note that the expression of Leibniz's rule (|3.20|) is exactly the same as for discrete deriva- 
tives j^]. For the sake of brevity we will write D x f := [A x ,f]. The expression f|3.20|) is 
operatorial. If we want to have a functional expression we need to project it by acting on 
a constant function as 1. In this case, however, the Leibniz rule ()3.2()j) is an identity. 

Taking into account the Leibniz rule ()3.20j) we can prove for the three g-derivatives 
introduced above the following property 

[A x , ([3 x x) n ] = D x ((3 x x) n = n{f3 x x) n -\ e N, (3.21) 

thus showing that {(3 x x) n are basic polynomials for the operator D x , and, when projected, 
for A x . So, we have defined an operator D x which on functions of j3 x x have the same 
properties as the normal derivatives d x on functions of x. Hence, we can say that what- 
soever is valid for differential equations can be also valid for the D x operators provided 
we consider instead of functions of x functions depending on f3 x x. This is the content 
of the q-umbral correspondence. In the case of linear differential equations, when the 
derivations act linearly on functions, the projection procedure will transform the operator 
D x into A x and we will get g-difference equations. 

Let us analyze the meaning of the basic polynomial operators (f3 x x) n for the three 
g-derivatives operators. Since 

(T x - l)x = x{q x T x - 1), (T x - \y l x = x{q x T x - (3.22) 

(1 - T- 1 )* = x(l - -T- 1 ), (1 - T~ l )~ l x = x{\ - -T- 1 )- 1 , (3.23) 

(T x - T x l )x = x(q x T x - -T-h), (T x - T" 1 )^ = x(q x T x - -T~ l y\(?>M) 

q x q x 



we get 



(3 x x = q x x(l + xd x )(q x T x - 1 
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p-x = q-x{\ + xd x )(l - -T- 1 )- 1 , (3.25) 

(? x x = q B x x(l + xd x )(q x T x - -T" 1 )- 1 . 

These expressions have an operational character. In order to reduce them to a function we 
have to project them by acting (3 x x on a constant. So, taking into account that T x 1 = 1 
and T~ l 1 = 1, in all the three cases we have 

P x x 1 = x. (3.26) 

One can demonstrate by induction that 

n ' 

(f3 x x) n 1 = (f3 x x)((3 x x) ■ ■ ■ (p x x) 1 = — -x n , Vn e N+, (3.27) 

[n] q \ 

where 

a n — 1 1 — a~ n 

W = K" = y^, (3-28) 



,2n 

I" 



q q zn - 1 1 

In'" 



, ( 3 - 29 ) 

q z — 1 g n g n 1 ^— ^ 



and 

[n] s ! = [n]„[n - • - • [1]„. (3.30) 
Consequently, if we consider an entire function like the exponential, we have 



00 \n 00 n 

Therefore, by the g-umbral correspondence, the exponential function becomes 

e A. = y A n^ _^ e A& n = y- A « * (3.32) 
^ n! ^ \n]J 

Let us consider the gaussian function, that we will be using later in this work. It takes 
the form 



oo 



-X(PxX) 



e — i = J2 ^^W"i = E ^r^fwrr (3 - 33) 

^— ' n\ n\ \2n\J 

By considering an arbitrary point xo and an arbitrary real constant a we have 

D x (p x x + x ) a = a(P x x + x ) a - 1 . (3.34) 
The proof of eq. ()3.34|) is based on the idea that the differential equation 

(x + x )d x f = af, (3.35) 



2n 
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whose solution is f — (x + xo) a , can be transformed by the g-umbral correspondence into 
the discrete equation 

((3 x x + x )D x (/3 x x + x ) a = a(f3 x x + x ) a , (3.36) 

which has formally the same solution in power series (substituting x by (3 x x in the discrete 
case). Effectively, 

00 a— nrrn— 1/ ,\ 

( x + XQ y = y x o U*^ fc V (3.37) 

n=0 

while 

(P x x + x o y = V X ° J (^) w (3-38) 

After projection we get 

(fcx + x o y 1 = jr Xl ~ nm =f a ~ k \ p x xT 1 = £ ^IB^^. (3.39) 

n=0 n=0 L iq 

□ 



4 Umbral functions 

In this section we consider some basic discrete functions obtained by means of the q- 
umbral method and discuss their range of validity taking as reference the corresponding 
solution of the difference equation and the continuous corresponding function. As concrete 
examples we shall consider as g-derivative operator two of the g-derivatives of the previous 
sections: the right A + and the symmetric A s . We present here two of the functions which 
appear in the symmetry reduction of the heat equation, the exponential and gaussian 
functions (J3.32j) . (J3.33|) . which exemplify the kind of results one can obtain from the q- 
umbral calculus for entire functions. These functions are described, respectively, by the 
ordinary differential equations: 

dfe(x) 

dx 

and 

dfg(x) _ 
dx 

Eq. (j4.1|) has as solution the function f e (x) = ae Xx and eq. ()4.2|) has as solution the function 
f g (x) = be~ Xx , where a and b are arbitrary integration constants. 



= A/ e (x) (4.1) 
-2Xxf g (x). (4.2) 



4.1 q— Exponential functions 

From the g-umbral correspondence the difference equation satisfied by the g-exponential 
is (see eq. (j4.1|) ) 

AE g (x) = XE q (x). (4.3) 
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We will discuss always the domain x > for all real A. The negative values x < can be 
obtained by changing the sign of A. 

The right exponential 

From (|3.1J) the difference equation (|4.3J) becomes 

E g (qx) = [l + {q-l)Xx]E q (x). (4.4) 

The solution can be expressed clS cL product 

n-l 

£,(? n *o) = II + (9 - 1) V'^o] ^o), n £ N+ (4.5) 

3=0 

We have four different behaviours of the g-exponential function according to the values 
of q and A, namely q j> 1 and A ^ 0. 



g > 1 



Let us at first consider the case A > 0. The recurrence of the difference equation 
implies that E q (x) is an (monotonous) increasing function in x. 

By the g-umbral correspondence the solution of the equation (j4.1j) can be also 
obtained by g-umbralizing the series representation of the exponential function. 
In such a case we have 



= E u [q " • (4-7) 



This solution converges for all x > 0, so that it gives the unique solution 
(l4.5H4.h1) of UQJ) in all the domain. 

— In the case A < we see that the recurrence leads to a decreasing function 
as long as 1 — (g — l)Xx > 0. However, for further values of x, where 1 — 
(q — l)Xx < 0, the function oscillates with higher diverging amplitudes. For 
the particular point x = l/((q — 1)|A|) we have E g (qx ) = 0, therefore also 
E q (q n Xo) = 0, Wn £ N. This means that in the g-lattice x n = xoq n ,n £ Z, 
the g-exponential decreases for n negative, i.e. for values of x less than xq and 
vanishes after, for n positive, avoiding the oscillations. 

As q — > 1 the point Xq — > oo, so that the g-exponential in that lattice becomes closer 
and closer to the (continuous) exponential. The g-umbral function (|4.7|) displays 
these same features as it is shown in Fig. 1. 
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1 2 3 4 5 



Fig. 1A. Plotting of the exponential function (continuous line) and right q- 
exponential function (dashed line) for A = 1 and q = 1.3. 




Fig. IB. Plotting of the exponential function (continuous line) and right q- 
exponential function (dashed line) for A = — 1 and q = 1.3. 
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Fig. 1C. Enlarged plotting of the exponential function (continuous line) and 
right q-exponential function (dashed line) for A = — 1 and q = 1.3. 

• q < 1 

— For A > and for big values of n one can always find [1 + (q — l)\q n xo] > and 
the solution is a monotonous increasing function for small values of xq. For 
the value xq such that [1 + (q — l)Axo] = 0, the solution is not defined, in fact 
it diverges in that point, and, consequently, also diverges in the lattice Xoq~ n , 
n G N. For the rest of the points x = xoqi such that [l + (q — l)Ag J a;o] < 0, that 
is, for the values x > xq , the solution oscillates, changing the sign alternatively, 
with divergences. The amplitude of the oscillations tends to zero, as x — > +oo. 

— For A < the solution is a monotonous decreasing function that in the limit 
x — > goes to zero. 

For the g-umbral series the radius of convergence is given by R — (\q — 1| |A|) -1 . 
The g-umbral function E q {x) has two vertical asintotes at the symmetric points 
x = ±R. Therefore, inside the range of convergence (where there are no 
oscillations) the g-umbral function reproduces correctly the solution E q (x), 
but for x > \R\ where it does not converges it supplies no information about 
the solution. These features are illustrated in Fig. 2. 




I , , , , , : r- . , , , , A 

2 4 6 8 

Fig. 2A. Plotting of the exponential function (continuous line) and right q- 
exponential function (dashed line) and the solution of eq. I|4.4j) given by the 
points, for A = — 1 and q = 0.5. 
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Fig. 2B. Plotting of the exponential function (continuous line) and right q- 
exponential function (dashed line) and the solution of eq. (|4.4|) given by the 
points, for A = 1 and q = 0.5. 
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Fig. 2C. Enlarged plotting of the solution of eq. (|4.4j) . for A = 1 and q = 0.5. 

The symmetric exponential 

When we consider the A s operator (jH.Hj) . the difference equation (|4.H|) becomes a three 
term relation given by 

E s q (qx) = {q- l/q)\xE s q (x) + E s g (x/q). (4.8) 

So, there are two independent solutions for the symmetric g-exponential. It is not neces- 
sary to distinguish the cases q > 1 and q < 1 because they play a symmetric role. For 
A > relation ()4.8|) gives growing functions in the x variable. For A < the function 
initially is decreasing but after a certain point the discrete solutions start to oscillate. 

The g-umbral solution to the recurrence equation ()4.8|) is given by 



w = y {x f iq ~ 1/q)k 



(4.9) 
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This series converges for all x, so that it provides one of the solutions to the recurrence 
([4.8)1 . and in the limit q — > 1 it goes into the continuous exponential. The behaviour of 
the umbral solution E* presents the features of the recurrence above described, i. e., it 
approaches the continuous exponential up to the first zero (for A < 0) but it becomes 
wildly oscillating beyond that point (see Fig. 3). 




Fig. 3A. Plotting of the exponential function (continuous line) and symmetric 
g-exponential function (dashed line), for A = 1 and q = 1.3. 



1 




1 2 3 4 5 6 



Fig. 3B. Plotting of the exponential function (continuous line) and symmetric 
g-exponential function (dashed line), for A = —1 and q = 1.3. 



14 



.001 
. 00075 
.0005 
. 00025 



-0 . 00025 
-0 .0005 
-0 . 00075 
-0 .001 




Fig. 3C. Enlarged plotting of the exponential function (continuous line) and 
symmetric g-exponential function (dashed line), for A = — 1 and q = 1.3. 

Comparing Fig. 1 and Fig. 3 we can see that the symmetric exponential function gives 
a better approximation than the right exponential. However, it can be easily shown that 
there is no initial condition for equation (|4.8j) such that the symmetric g-exponential 
vanishes for all subsequent points, as it was the case for the right exponential. So, there 
is no way to avoid the oscillations. 



4.2 q- Gauss ians 

In this case eq. ()4.2|) . which has as a solution the gaussian function, becomes the in- 
difference equation 

AG q , x (x) = -\xi3 x G q>x {x). (4.10) 

In the following we will discuss briefly the cases of the right and symmetric g-gaussians. 
The right gaussian 

The difference equation (|4.1()jl becomes a three-term recurrence equation 

cT'GUfx) - (q- 1 + l)G+ >x (qx) + G+ A (x) = -X(q - l) 2 x{x 2 d x + 1)G+ A (x). (4.11) 

However, now eq. (|4.11|) is a differential-difference equation. So, it is quite difficult to 
find directly the solutions or even to discuss the general behaviour of its solutions. We will 
assume, as shown in the case of the exponential function, that, whenever the g-umbral 
series is convergent it will converge to the solution of the difference equation. 

The g-umbral series supplies a solution to eq. (|4.11|) . given by 

^^HfWh^. ( , 12) 

For q > 1 the series converges for all x, but for q < 1 the series diverges everywhere (for 
i / 0). Some plottings of G^ x are shown in Fig. 4. The behaviour is similar to that 
of the g-exponential. Whenever we have a decreasing function of x, at a certain point 
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it vanishes and beyond that point it starts to oscillate with increasing amplitudes, thus 
departing from the behaviour of the continuous gaussian function. 
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Fig. 4A. Plotting of the gaussian function (continuous line) and right q- 
gaussian function (dashed line), for A = 1 and q = 1.3. 
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Fig. 4B. Enlarged plotting of the gaussian function (continuous line) and right 
g-gaussian function (dashed line), for A = 1 and q = 1.3. 
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The symmetric gaussian 

The recurrence relation of eq. (|4.1()jl for this case is an even more involved differential 
difference equation than eq. (|4.11|) . so we prefer not to write it down. The g-umbral series 
solution is given by 

-\x 2 ) k (2k)\(q- l/q) 2k 



G s q , x (x) 



E 



The radius of convergence is R = oo. For the symmetric gaussian function we have similar 
results as for the right gaussian, which, however, are valid for any value of q (see Fig. 5). 



(4.13) 
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Fig. 5A. Plotting of the gaussian function (continuous line) and symmetric 
g-gaussian function (dashed line), for A = 1 and q = 1.3. 
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Fig. 5B. Enlarged plotting of the gaussian function (continuous line) and 
symmetric g-gaussian function (dashed line), for A = 1 and q = 1.3. 
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In conclusion, we can say that the decreasing asymptotic behaviour of the classical 
functions is not fully reproduced by the corresponding umbral g-functions. In that region 
the discrete functions approach the classical ones up to a point where they vanish, but 
beyond this point they oscillate going far away from the continuous analogues. Therefore, 
a good parameter measuring the radius of the domain where the g-functions imitate the 
continuous functions is given by the first zero in the region of asymptotic behaviour. This 
is depicted in Fig. 6 for the g-exponentials and the g-gaussians. 

From these plottings we evince that the domain of convergence of the g-exponential 
function to the exponential function increases in a monotonic continuous way as g — > 1 
for the right exponential while in the symmetric case there are discontinuities for small 
values of q. The situation is slightly different in the case of the gaussian function, when, 
as q decreases also the domain decreases up to a minimum A dependent value go- Below 
go the domain increases as q — > 1. In the case of the symmetric gaussian function we 
again have that for q small the function is discontinuous. 
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Fig. 6A. Plotting of the position of first zero of the right g-exponential as 
function of q for A = — 1 
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Fig. 6B. Plotting of the position of first zero of the symmetric g-exponential 
as function of q for A = — 1 
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Fig. 6C. Plotting of the position of first zero of the right g-gaussian as function 
of q for A = 1 
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Fig. 6D. Plotting of the position of first zero of the symmetric g-gaussian as 
function of q for A = 1 

5 Examples 

5.1 Discrete heat q— equation 

Taking into account the g-umbral correspondence 

d x -> A x , x -> (3 x x (5.1) 

we obtain from any linear differential equation with constant coefficients an operator 
equation which, when projected, gives us a g-discrete equation. In the case of the heat 
equation we get 

(d t -d 2 xx )u = =}► {A t -A xx )u = (5.2) 

Let us consider now the problem of obtaining the symmetries for the g-discrete eq. (|5.2|) . 
We can apply the g-umbral correspondence also to the determining equations, as they 
are linear in the infinitesimal coefficients £, r and /. So, making use of the Leibniz rule 
we obtain the following set of determining equations: 

D x {t) = 0, (5.3) 

A(r) - 2Q*(£) = 0, (5.4) 

- Q xx {® - 2W) = 0, (5.5) 

D t (f) ~ Q xx (f) = 0, (5.6) 

where D x {r)l = [A x ,r]l = A^r. Moreover 

= D x (D x (f))l = [A x , [A X J]]1 = A X AJ. 

Note that these determining equations have formally the same expression for all q- 
derivatives, for the continuous derivatives and in the discrete case studied in From 
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eq. (|5.ip . the solution of this system ()5.3|) is 

r = r 2 ({3 t t) 2 + n(f3 t t) + r , (5.7) 
e = \{n + 2T 2 (p t t))(p x x) + U&t) + £„, (5.8) 

/ = ^(/U) 3 + \r 2 {(3 t t) + + 7, (5-9) 

where r , r 1; r 2 , £o ; £i and 7 are arbitrary functions of T x , T t and of q x and g t . 

By a suitable choice of the functions n, r 2 , r , £i, £o and 7 we get the following 
representation of the symmetries: 





(A t u)d u , 


(5.10) 


p? = 


(A x u)d u , 


(5.H) 


w q = 


-- ud u , 


(5.12) 


B q = 


(2(/3 t t)A t u + ({3 x x)A x u)d u , 


(5.13) 


D q = 


(2(p t t)A t u + {(3 x x)A x u + -u)d u , 


(5.14) 


K q = 


{{(3 t t) 2 A t u + {p t t){{3 x x)A x u + -^ x x) 2 u + -(p t t)u)d u , 


(5.15) 



that close into a 6-dimensional Lie algebra, isomorphic to the symmetry algebra of the 
continuous heat equation. Another realization of this algebra was obtained in ^H], by a 
different procedure and used to find symmetric solutions of the discrete heat equation. 

Now we use these symmetries to construct few solutions of the g-discrete equation 
(J5.2|) . Taking into account the symmetries Pq and P q , choosing the g-parameters in such 
a way that q x — q t — q and, using the variable separation method [20], we can write the 
solution of the g-heat equations as 

u(ptt,p x x)=v(ptt)w(p x x). (5.16) 

Then eq. ()5.2|) reads 

(A t v(p t t))w{p x x) - v(f3 t t) (A xx )w(/3 x x)) = 0. (5.17) 
From eq. ([5.17|) we deduce that the functions v and w must satisfy the following equations 
A t v(p t t) = \v(p t t), A x w(/3 x x) = \f\w(p x x). (5.18) 
By means of eq. ()4.3|) . we have 

v(t) = v(/3 t t)l = e^l = riT = < ( 5 - 19 ) 

n=0 [ iq - 
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and 



X n/2 



w(x) = w(p x x)i = e ^ x i = = e f x - ( 5 - 2 °) 

Hence, the solution will be 



n=0 1 iq 



u(t,x) = e*ef x . (5.21) 

Let us consider now the symmetry reduction with respect to the operator B q ^H]- In 
this case introducing the appropriate symmetry variable 77 = -%= we get 

V ftt " 

<^m^ Ji m^ (5 - 22) 

The solution (|5.22|) of the g-heat equation is meaningful as long as we are considering 
positives times and the value t = is out of our time domain. In such a situation the 
solution (|5.22jl is entire and can be represented as a Taylor series and, thus, g-functions 
like the gaussian or the square root are meaningful. The method would not provide a 
meaningful g-function if we would consider all values of t. However, in t = also the 
boost solution of the continuous heat equation would be singular and, thus, meaningless. 



5.2 A generalized Hermite equation 

Let us consider the following g-difference equation 

A xx ^j + (3 x xA x ^j = E^j. (5.23) 

By explicitating the operator f3 x appearing in eq. ([5.23p according to (|3.11|) . we can rewrite 
it as 

A xx tP + xtp x = Eip. (5.24) 

By the g-umbral correspondence eq. ()5.23J) is related to the ordinary differential equation 
ipxx + xipx — Eip whose solutions are given in terms of gaussian functions and Kummer 
confluent hypergeometric functions. Explicitly we have [2*T] : 

tP{x) = A ie -^M(^ + ^E,^x 2 )+A 2 xe-^M(l + ^E,~x 2 ) (5.25) 

= e 2 <p{x). 

Both the gaussian and the Kummer confluent hypergeometric functions have a power 
series expansions in terms of x, 

4>{x) =A l [l + (l+ \e)^ + (1 + \E){?> + \E)^ + ...]+ (5.26) 

A 2 [x +(2 + -E)- + (2 + -E)(A + -E)- + . . .]. 

So, the g-differential difference equation ()5.24|) will have a solution in terms of g-gaussian 
and g-Kummer confluent hypergeometric functions. Let us notice that, if E = —2n with 
n a positive integer number, the g-Kummer confluent hypergeometric function is just a 
polynomial. 
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6 Conclusions 



In this paper we present a g-extension of the umbral calculus and use it to provide 
solutions of linear g-difference and g-differential difference equations. In this way we 
obtain solutions which have a continuous limit. 

The discretization procedure given by the recipe d x — ► and x — ► j3 x x works well 
for linear equations also in the case of g-shifts operators. In particular, it preserves the 
classical Lie symmetries which are described by linear equations. 

We study in detail the behaviour of the g-exponential and g-gaussian functions and 
show their range of validity which depends from the g-discrete delta operator under 
consideration. The domain of convergence of the g-function to the continuous function 
is characterized in term of the zeroes of the g-function. The results are usually better in 
the case of symmetric g-delta operators. 

Further work is in process on the complete description of a coherent g-umbral calculus 
and a comparison of the discrete and g-discrete solutions. 
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